Abstract. We present times series of high resolution 
spectra of AGB variables at 4 /jm. Line profiles from the 
major contributors to the spectra of oxygen rich stars at 
4 /mi, OH, H2O, HC1 and SiO, are examined. The veloc- 
ity as well as shape variations of these profiles with time 
are discussed. The line profiles investigated frequently 
have emission and multiple absorption components. The 
changes with time of the 4 /xm region lines do not always 
follow the cyclic variability seen in NIR spectra and in 
the photometric light curve. We interpret and discuss the 
y—{ . results qualitatively considering comparing the spectral 
^^ ' variability with that of the well behaved 1.6 /xm region 
^^ . and of dynamical model atmospheres. Miras and semireg- 
ular variables are compared. The origins of non-periodic 
behavior are discussed, including the role of spatial inho- 
mogeneities in the stellar atmosphere. 
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1. Introduction 

The structure of the envelopes of Asymptotic Giant 
Branch (AGB) stars is dominated by large amplitude stel- 
lar pulsations. The mechanical energy of the pulsation 
can b e larg er than L/c by orders of magnitude (Hinkle 
et al. 1982| ), resulting in a cool, highly extended outer 
atmosphere. As a result of the pulsational flow of me- 
chanical energy, the structure of the atmosphere varies 
strongly during the light cycle. It is possible with time se- 
ries spectra to monitor the flow of material in the stellar 
atmosphere. In addition, velocities measured from differ- 
ent spectral transitions of different atomic and molecular 
species probe the dynamics at different (optical) depths 
in the atmosphere. 

One result of the extended, cool atmosphere of these 
stars is that the stellar flux is primarily emitted in the 
infrared. Due to the rich infrared molecular spectra of the 
cool AGB stars, the infrared is also the optimum spectral 
region for studying the outer atmosphere of these stars. 
The analysis of time series of high resolution infrared spec- 
tra has been used over thre e decades b y Hink le and collab- 
orators (e.g. Hinkle et al. |1982[ |1984| , |1997| ). Both miras 
and the lower amplitude semiregular variables have been 
investigated. The authors mainly used CO lines in the 1.6- 
2.5 /im region, since those transitions can be easily identi- 
fied and observed. The high excitation first overtone CO 
lines and the second overtone lines of CO show velocity 
variations of 25 to 30kms _1 in miras and about lOkms -1 
and less in semiregular variables. In the 1.5-2.5 /im region 
other atomic and molecular lines have been found to have 
similar behavior, while low excitation first overtone lines 
of CO show only very small velocity changes with phase. 
It has been suggested that these spectral lines as well as 
some of the H 2 lines probe a more extended pseudo- 
stationary layer. Spectral lines with excitation potentials 
near zero volts and large column densities, e.g. the lowest 
excitation CO fundamental lines situated around 4.6 /im, 



probe the classic, expanding circumstellar shell (Ryde et 
al. p000| ). 

In this paper we will investigate the molecular spec- 
tra in the 4/im region. In oxygen rich atmospheres, the 
principal molecular constituents are the first overtone SiO 
bands, the fundamental bands of OH, and the HC1 funda- 
mental. In addition, the 4 /xm region is blanketed by a large 
number of weak water lines in stars with 0>>C. We will 
discuss high spectral resolution observations of lines from 
all four of these molecules in this paper. This is the first 
use of these 4 /im lines at high resolution for monitoring 
the atmospheric changes in pulsating stars. However, the 
4 /im spectrum of SiO has been monitore d pre viously at 
medium and low resolution (Hinkle et al. 1976 ; Rinsland 
fc Wing |l982T) . 



Aringer et al. ( 1997 , 1999b ) have shown from medium 
resolution 4 /mi spectra that using hydrostatic models can- 
not reproduce the spectra properly. In order to match the 
observations one needs to introduce additional features 
like the "warm molecular envelope" proposed by Tsuji et 
al. ( 1997 ). On the other hand, calculations based on re- 
cently developed dynamical atmospheric models (Hofncr 
et al. 1998) ) have the ability to reproduce the observed 
spectral features as shown by Aringer et al. (1999b) for 
the SiO first overtone band heads around 4/im. In this 
paper we will present the observational data for a sample 
of miras and semiregular variables. These data will pro- 
vide the input for a forthcoming second paper (Aringer et 
al. in preparation) that will deal with a comparison of our 
findings with atmospheric models focusing on the behav- 
ior of the SiO lines as in the previous papers by Aringer 
et al. (|1997| 7 |l999bt). 
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2. Molecular data 

2.1. OH 

The hydroxyl free radical is a well known contributor of 
vibration-rotation lines to the spectra of K and M stars. 
Aringer ( pOOO ) finds that OH appears in giant spectra 
at an effective temperature below approximately 4200 K. 
The OH molecule is well known for the three maser lines 
in the rotational spectra of AGB stars at 1612, 1665 and 
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1667 MHz, respectively. These lines provide an important 
tool for deriving accurate center-of-mass velocities of AGB 



the band heads at medium spectral resolution, Hinkle ct 



stars (see Habing 1996 for a review). Within our sample 
of stars OH masers have been detected only in R Cas and 



RLeo (Benson et al. 1990). The maser lines are much lower 
energy transitions than any of the OH lines discussed in 
this paper and are formed in the expanding circumstcllar 
shell. 

The 3-4 (im region contains the vibration-rotation fun- 
damental bands of OH. The first overtone vibration- 
rotation bands are found around 1.5 fim. Due to the low 
molecular weight of OH, the rotational constant is large 
and OH does not form band heads in the P branches in 
the near infrared. The OH lines are out of the X 2 n ground 
state and exhibit large lambda doubling. Hence the lines 
appear in quartets in the spectrum. A detailed study of 
the first overtone OH lines in the mira variable RLeo is 



described by Hinkle (1978). Line positions derived from 
Coxon (1980) and Coxon & Foster (1982) are used in the 



current paper. The line list consists of a group of 17 1-0 
and 2-1 P branch lines selected from the total of 46 lines 
in the 2500 to 2800 cm -1 region. The chosen lines have 
minimal blending with telluric and other stellar features. 
The entire set of 46 lines was used to evaluate the mea- 
surement process as described below. 



2.2. SiO 

SiO is, after CO and H2, one of the most abundant 
molecules in the atmospheres of oxygen- rich AGB stars. It 
is formed in all cooler oxygen rich stars. Model spectra in- 
dicate that significant amounts of SiO are formed in stellar 
atmospheres with effect ive te mperatures of approximately 
4000 K or less (Aringer 2000 ). In addition to being abun- 
dant, SiO is a primary constituent of oxygen-rich dust. 
Since dust is one of the requirements for AGB mass loss 
(Gail & Sedlmayr 1986 ), SiO formation is a prerequisite 
for mass loss. 

SiO maser emission is very common among miras (Jew- 
ell et al. 1991). Beside sharp maser lines, broad thermal 
SiO emission in rotational lines has been detected in sev- 



eral AGB stars (see Habing 1996 for a review). Some SiO 



maser lines involve relatively large J rotational transitions. 
These lines could be formed near the same region probed 
by the infrared vibration-rotation SiO lines. Two infrared 
SiO vibration-rotation transitions have been observed in 
late-type stars. The fundamental bands are found in the 
spectrum near 8 /mi. The first overtone bands are found 
near 4 /im. We report here on measurements of lines from 
the first overtone. 

Unlike OH, the SiO line spacing is small. This spacing 
combined with the transitions occurring in a X 1 !] state re- 
sults in well defined R-branch band heads. The first high 
resolution spectra of the SiO fi rst ov ertone lines in a star 
were obtained by Beer et al. ( |l974|) who used the spec- 
tra to derive Si isotopic abundances in a Ori. Observing 



al. (1976) found that long period variables show a strong 
variability in the intensity of the SiO bands at 4 /mi. In 
miras, around light maximum SiO becomes almost unde- 
tectable, strengthening towards light minimum. This re- 
sult has been confirmed by Rinsland & Wing ( 1982Q , who 
also noted that SiO absorption features are most promi- 
nent in stars with the lowest temperatures (i.e. at light 
minimum in miras). Different explanations for the strange 
behavior of SiO in AGB variables have been discussed in- 
cluding dissociation, additional absorpti on (e .g. by H2O) 
or emissio n of SiO (Rinsland & Wing 1982 , Aringer et 
al. 1999b). The only published high resolution spectra of 



the first overtone bands of SiO in AGB variab les is an 
atlas of the 4 /mi range by Ridgway et al. ( 1984 ). 



The SiO line list used in this paper is made up of 20 
unblended 28 SiO lines in the 2-0 R branch between the 
3-1 band head and the 2-0 band head (2-0 R22 to R46) 
spanning the wavenumber interval 2475 to 2493 cm -1 . The 



frequencies used are those reported in Hinkle et al. (1995) 



2.3. HCl 

There are no previous detailed investigations of HCl in 
cool evolved stars. Ridgway et al. ( 1984 ) report lines of 
HCl in the 4 /im spectrum of the S-type mira R And, but 
they did not find HCl in oxygen- or carbon-rich AGB 
stars. The HCl vibration-rotation fundamental crosses the 
3-4 /im region. The transition is a X 1 !] and since HCl, 
like OH, is a hydride the rotational constant is large and 
hence the line spacing is large without prominent band 
heads. The frequencies of the lines are taken from Le 
Blanc et al. p9|). Lines of both the H 35 C1 and H 37 C1 



isotopes were used. This augments the rather spare num- 
ber of lines present and since the solar system 35 C1/ 37 C1 
is small (about three) the intensities of the lines from the 
two isotopes should be similar. The line list includes 20 
lines with the least blending selected from the 1-0 and 2-1 
bands limited to the wavenumber region between 2500 and 
2800 cm- 1 . 

2.1 H 2 

Water was first observed in the infrared spectra of cool 



stars in the early 1960s (see Spinrad & Wing |l969[ ). Later 
in that decade observations showed water bands to be 
very strong, variable features in mira spectra (Johnson 
& Mendez |l970p . Recently Alvarez et al. fl2000| ) con- 



cluded that the detection of deep water bands in stellar 
spectra is a direct indication of pulsation. Early work- 
ers also recognized that despite water's large abundance 
in cool, oxygen-rich stars, observations of this molecule 
from the ground are inhibited by water bands in the earth 
atmosphere. Water m ust either be o bserve d from space 
(Arin ger et al. 1999a Matsuur a et a l. 1999 ; Truong-Bach 
et al. 1999J ; Yamamura et al. 1999 ) or by using excited 
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levels which are populated in atmospheres of stars but 
are not populated in the Earth's atmosphere (Hinkle & 
Barnes |l979| ) . 

In addition to the vibration-rotation infrared spectrum 
of H 2 0, this molecule has a rich pure rotational spec- 
trum in the microwave and far-infrared. While consider- 
able blockage of this spectrum also occurs as a result of 
water in the Earth's atmosphere, water maser lines are 
known in circumstcllar envelopes of evolved stars (Mentcn 
& Melnick 1989) including three of our program stars, 
namely RCas (Colomer et al. |2000 ), RLeo (Menten k, 
Young |1995| ), and RXBoo (Colomer et al. |2000D . 

The spectrum of water in the atmospheres of LPVs 
originates in two distinct regions of the star's atmosphere: 
Most of the lines are formed in the circumstellar re- 
gion closest to the star, while some lines follow the pho- 
tosphcric behavior and therefore originate much deeper 



Table 1. Observing log for the miras. The period in col- 
umns is taken from the GCVS. 



in the atmosphere (Hinkle & Barnes 1979; Matsuura et 
al. 1999] ). Infrared H2O absorption cannot originate far 
out in the circumstellar envelope as this would not be 
consistent with the mass loss rates found in these objects. 
Tsuji et al. ( |1997| ) and Tsuji ( pOOOj ) argue that the non- 
photospheric water lines indicate the existence of a warm 
molecular envelope. 

In this paper we will discuss high excitation lines of the 
V3 vibration-rotation band which are present in the 4 fim 
spectrum. The band origin is at 3756 cm -1 . The line list 
consists of 102 high excitation 1/3 lines between 25 00 an d 
2800cm" 1 . The frequencies are from Zobov et al. ( 2000|) . 
The line list used for this paper involves moderately high 
rotational (J) quantum numbers between 12 and 24 and 
should not include lines formed in the inner circumstcl- 
lar region. However, we note that like the OH and HC1 
transitions, the 1/3 transition is a fundamental vibrational 
transition. Fundamental transitions are the most likely to 
be detected in cool conditions and are also the preferred 
routes for radiation resulting from recombination. 

3. Observations 

The spectra were obtained by one of us (K.H.H.) with 
the Fourier Transform Spectrometer (FTS) operated at 
the Mayall 4m telescope at Kitt Peak National Obser- 
vatory. A description of the instrument can be found in 
Hall et al. (1979). The FTS is no longer available for use, 
having been closed in the mid-1990's. The spectra are 
an inhomogeneous set, having been observed for several 
projects. However, all cover the 2400 to 2800 cm -1 region 
required by the line lists. The resolution achieved is about 
R=100000 except for very few observations at a somewhat 
lower resolution but none below 50000. The observations 
were obtained mainly between 1984 and 1986 with a few 
additional spectra observed earlier. Three spectra of x Cyg 
were observed in 1994 and 1995, shortly before the FTS 
on Kitt Peak was closed. All d ata have been reduced as 
described in Hinkle et al. (|l982|) . 
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JD 


Visual 


S/N 
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type 




2440000+ 


phase 




RCas 


Mira 


430 


5986 


0.54 


53.7 








6018 


0.61 


106.7 








6108 


0.82 


27.7 








6160 


0.94 


87.9 








6192 


0.01 


47.7 








6223 


0.09 


52.1 








6243 


0.13 


35.1 








6312 


0.29 


92.6 








6454 


0.62 


44.0 








6510 


0.75 


50.2 








6570 


0.89 


58.4 


Cet 


Mira 


332 


3437 


0.89 


243.4 








4775 


0.81 


108.6 








5986 


0.45 


63.2 1 








6018 


0.55 


93.4 








6095 


0.78 


2 








6108 


0.82 


57.6 








6312 


0.43 


66.6 


X Cyg 


Mira 


408 


3437 


0.17 


135.9 








5986 


0.35 


53.8 








6018 


0.43 


88.0 








6108 


0.65 


75.0 








6160 


0.77 


86.3 








6192 


0.85 


85.8 








6223 


0.93 


65.2 








6243 


0.98 


36.0 








6312 


0.15 


71.4 








6346 


0.23 


59.4 








6427 


0.43 


22.0 








6454 


0.50 


50.6 








6510 


0.63 


64.4 








6570 


0.78 


65.4 








9323 


0.62 


89.4 








9363 


0.71 


93.7 








9450 


0.93 


88.0 








9612 


0.30 


120.1 








9798 


0.75 


66.3 


R Leo 


Mira 


310 


4650 


0.60 


329.2 








5986 


0.76 


70.3 








6164 


0.31 


54.6 








6427 


0.15 


24.4 








6511 


0.41 


74.6 



For four miras the data set includes time series of 
4/im spectra, allowing the stellar pulsation to be moni- 
tored. There are also single epoch observations of a num- 
ber of semiregular variables. Together with some funda- 
mental parameters from the General Catalogue of Vari- 
able Stars (GCVS, Kholopov et al. |1985 
vational details are given in Tables 






the obser- 
Phases are 



and~p 
calculated from t he closes t light maxima provided by the 
AAVSO (Mattei |l984-95|) . As can be seen from Tabled 
for three of the four miras (RCas, oCet and xCyg) ob- 
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Table 2. Observing log for the semiregular variables. The 
period in column 3 is taken from the GCVS. Phases are 
rough estimates from AAVSO light curves. 



GCVS 


Var. 


Period 


JD 


Visual 


S/N 


name 


type 




2440000+ 


phase 




RZ Ari 


SRb 


30 


7071 


- 




RXBoo 


SRb 


340 


6986 


0.0 






WCyg 


SRb 


131 


6160 


0.4 




42.8 


RUCyg 


SRa 


233 


6160 


0.2 




16.9 


g Her 


SRb 


89 


6986 


0.9 




i 


RLyr 


SRb 


46 


6986 


- 




SW Vir 


SRb 


150 


6986 


- 





servations from different light cycles but at similar phases 
exist. This can be used to investigate the cycle to cycle 
behavior of these stars (see below) . The dust enshrouded 
object IRC+10216, which is a featureless, bright contin- 
uum source at 4/im, was also observed. The IRC+10216 
spectrum allows easy identification of telluric lines. 

Seven stars of our sample (all miras, RU Cyg, W Cyg, 
gHer) have previously been investigated using 1.6 /im 
spectra (Hinkle et al. |1982| , |1984| , |1997| and Lebzel- 



1999]). These papers also give a summary of the fun 



tcr 

damental properties for these stars from the literature. 
Three miras of our sample are 'typical' representatives 
for oxygen rich disk miras with respect to period, ampli- 
tude, spectral type and mass loss. Since they are among 
the brightest miras in the sky, these objects are very well 
studied. The fourth mira is x Cyg, a S-type star, the pul- 
sational beha vior o f which has been studi ed in detail by 
Hinkle et al. (|l982|) and Wallerstein (|l985| ). 

The spectra are very blended and it is typically not 
possible to find 'clean' individual line profiles. Hence we 
have worked from average line profiles calculated from all 
lines of the list. Velocities have been determined from the 
deepest point in the average profile of the lines for each of 
the four molecules investigated. In many cases this veloc- 
ity summarizes a blend so, in addition, the average profiles 
themselves are discussed below. The line lists typical con- 
tain about 20 lines for each molecule of similar intensity 
and excitation. For water, many more lines were available 
and as can be seen in the average profiles the 'noise' was 
significantly reduced. 

The average profile is modified by blending lines 
around each of the individual line positions. Given the 
high level of blending of some spectra it is remarkable 
that a reasonable average profile emerges. In this light it 
is critical that the average line profiles not be over inter- 
preted. We feel that the general features, especially those 
seen in time series and in the profiles of other molecules, 
are real. However, detailed features seen only in one pro- 
file could be artifacts of the averaging process. In the case 
of OH we performed a test where average profiles were 
made up from two independent line lists of 23 lines each. 



There was remarkable agreement between the two profiles, 
indicating that essentially all the features are real. 

The velocity of telluric lines was measured to check 
on the zero point of the velocity scale. This correction is 
typically done in FTS spectra because the reference and 
stellar beams can have non-parallel paths. In each of the 
4 fim spectra we found the telluric velocities to be dis- 
tributed around zero with a full width at half maximum 
of this distribution of about 0.1km s _1 . At longer wave- 
length alignment errors become less significant and this 
seems confirmed by this data set. Applying the small, ap- 
parently random telluric correction would increase the un- 
certainties and hence we have not done so. In any case the 
correction is less than or on the order of ±0.2 km s _1 in the 
measured velocities. Further details on 4m FTS data and 
FTS data reduction can be found in Hinkle et al. (|1984|, 



19951) 



4. Mira Results 

Time series of averaged SiO, H2O, HC1 and OH lines of 
RCas are plotted in Figs. |l] and ||, respectively. Typi- 
cal line profiles from the other three miras can be found 
in Figs. @ to pi No continuum is marked in the figures 
as the true continuum is at a much higher level due to 
the almost continuous absorption by weak H2O lines (and 
similar molecules) in this spectral region. The line pro- 
files shown are averaged line profiles so that the range to 
both sides of the line can be used to define some kind 
of pseudo continuum. In Figs. |l] to [| spectra have been 
normalized to this pseudo continuum. The center-of-mass 
velocity, for all objects from radio observations of thermal 
CO, is marked in the figures (references on the center-of- 
mass velocity as well as previously measured infrared and 
visual velocities can be found in Hinkle et al. 1982, 1984, 



1997). Lines could be detected at all phases but vary con- 



siderably in shape and strength. A quick comparison of 
similar phases shows sometimes large variation in the line 
profiles between the same phase in different light cycles 
for the same star (e.g. compare the spectra obtained at 
phase 0.61 and 0.62 in Fig.nh. In the 4/im region RCas 
is an extreme example of cycle-to-cycle variations. The 
other miras of our sample show some differences between 
line profiles obtained from different light cycles at simi- 
lar phases but these differences are typically not as large 
as in RCas. Note that observations obtained at similar 
phase but several years apart can also have very similar 
profiles (e.g. Fig. g). This is suggestive of a strong aperiod- 
icity as well as a periodic variation in the spectra. This is 
a common comment made about visual spectra and some 
features (e.g. low excitation 2-0 CO lines) in 2 /im spectra 
but is not seen in the highly periodic behavior of the great 
majority of lines in the 1.6-2.5 /im region (e.g. Hinkle et 
al. |l982f ). 

Emission components are clearly visible in several ob- 
servations. This suggests a filling of the absorption with 
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40 

velocity [km/s] 

Fig. 1. Averaged SiO 2-0 (left) and H 2 r/ 3 (right) line 
profiles of R Cas. Date of observation and according phase 
are given in the plot. The dotted line marks the center-of- 
mass velocity. 



40 

velocity [km/s] 

Fig. 2. Averaged HC1 fundamental 1-0, 2-1 (left) and OH 
fundamental 1-0, 2-1 (right) line profiles of R Cas. Date 
of observation and according phase are given in the plot. 
The dotted line marks the center-of-mass velocity. 



emission as an explanation for the weakening of the SiO 
lines observed in low resolution spectra as proposed by 



Aringer et al. (1999b). Multiple absorption components 
with varying strength are an obvious characteristic of the 
investigated lines in most of the spectra. Unresolved mul- 
tiple emission components also appear. 

For a given observation the four molecules investi- 
gated, OH, H 2 0, SiO and HC1, differ quite strongly in 
the shape of their line profiles. OH, as seen in the RCas 
spectra, can show an emission line while the other three 
molecules show an absorption line. Among the studied 
oxygen rich miras differences of this type are largest for 
RCas and smallest for RLeo. In x^yg OH and H 2 
behave much more similarly than in the other stars of 
our sample as illustrated in Fig. pi Furthermore, the four 
molecules also differ in their variability. In the following 
we will now discuss the line profile variations for each 
molecule separately. We will offer an interpretation of the 



complex line structures seen keeping in mind the limita- 
tions of doing this without detailed model calculations. 

If not explicitly stated velocities and velocity shifts 
mentioned in the following are relative to the center-of- 
mass velocity of the star. The convention is used with 
negative velocity indicating outflow of material. The line 
profile figures have heliocentric velocity for the abscissa. 

4.1. OH 

Fig. pi gives the resulting velocities derived from OH lines 
for RCas and \ Cyg plotted against time. Time is selected 
for the abscissa rather than phase since the variability 
does not always follow the photometric period. Velocities 
were determined only for those absorption components of 
the line which could be clearly identified and attributed 
to the line investigated. As can be seen in Figs, [j] and || 
other weak absorption components can exist in the line 
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20 -20 

velocity [km/s] 

Fig. 3. Averaged OH fundamental (1-0, 2-1) (left) and 
H2O vj, (right) line profiles of xCyg. Date of observation 
and according phase are given in the plot. The dotted line 
marks the center-of-mass velocity. 




20 -20 

velocity [km/s] 





1 1 ' : 


1 1 




0.65 


/Y"\i 


Cet 


- 


0.60 






- 


0.55 


X__ — /-'' \ 




'- 


0.50 




\ / / 


- 


0.45 






- 


0.40 




\jl 


- 


0.35 


1 , : 


1,1, 


- 



velocity [km/s] 

Fig. 5. Averaged SiO line profiles from two observations 
of o Cet at similar phase but widely separated in time. The 
solid line is the profile from Feb 11 1985 ($ = 0.82), while 
the dash-dotted line marks the profile at Oct 21 1977 (<J> = 
0.89). The dotted vertical line indicates the center-of-mass 
velocity. 




velocity [km/s] 

Fig. 6. Same as Figureg for OH. 

profiles. Obviously for complex line profiles such as these 
the velocities are subject to interpretation. For instance, 
doubled absorption can result from an emission line filling 
the center of a broad absorption line. We base the results 
extracted from the averaged line profiles on similarities 
found for CO lines of different excitation published before 



Fig. 4. SiO & H 2 line profiles of R Leo. Same as Figure|l|. 



(e.g. Hinkle ct al. 1982). For \ Cyg only the measurements 
between 1984 and 1986 have been plotted to allow for the 
same time interval as RCas. 

For both stars we see a velocity curve resembling the 
velocity curves found for CO first overtone high excita- 
tion lines in the 2.3 /zm region (see e.g. Hinkle et al. 1982). 
However, a closer inspection reveals significant differences 
between the two stars and relative to the CO velocity 
curve. For RCas we have a phase of line tripling between 
JD 2446190 and JD 2446300 (see the individual line pro- 
files in Fig.0). The first part of the velocity curve shows 
a layer at rather constant infall velocity and a variable 
component changing from outflow to infall. During this 
time of infall, absorption is becoming visible at an out- 
flow velocity again. Surprisingly, the temporal separation 
between observations of similar velocity is shorter than the 
pulsation period of 430 days, in the one cycle examined 
it is ~300 to 350 days. However, this behavior is clearly 
not related to the general pulsation of the star since the 
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Fig. 7. Velocities measured from absorption components 

of OH lines in xCyg (upper panel) and RCas (lower 

panel). The dotted line indicates the center-of-mass ve- 
locity. 



AAVSO light curve of R Cas at the time of the infrared 
spectroscopy shows a cycle length of the typical 430 day 
photometric period. The only notable deviation from the 
normal light cycle is a very bright maximum preceding the 
time these spectra were observed. The first light maximum 
we monitored has been a slightly weaker than usual. The 
absorption at constant infall velocity is not visible in the 
second half of the studied time interval and there is also 
no indication of a further line doubling at the end of the 
second cycle around JD 2446600. 

In x Cyg OH shows velocity changes with the period 
of the photometric variations, i.e. observations of similar 
velocity are separated by roughly one period of the star. 
Between the two cycles we have one observation with ab- 
sorption line doubling. We also found line doubling on 
JD 2446749, but with a much smaller velocity difference. 
In addition, we note that the averaged line profiles of 
X Cyg indicate further phases of line doubling, but in these 
cases the second component is very weak not allowing an 
accurate velocity determination. 

Beside the absorption components OH emission is visi- 
ble at several phases and can be quite strong. In the R Cas 
spectra an OH emission component can be detected at all 
phases, but in several spectra it is not obvious since the 
emission peak is embedded in the absorption profile re- 
sulting in a weaker emission line that does not exceed the 
level of the continuum. The emission component starts 
at the blue end of the line at a velocity relative to the 



center of mass of roughly -14.5 km s _1 (i.e. +2kms~ 1 he- 
liocentric) and moves redwards. Around JD 2446312 the 
emission component reaches the center-of-mass velocity. 
Four months later a new emission line is visible at the 
blue edge of the line again. It is not clear whether the first 
emission line is still visible - now at an infall velocity - 
on JD 2446510. On the other hand, the second emission 
line can already be visible when the first one reaches the 
center-of-mass velocity. Unfortunately, the coverage of the 
variability is not very good during this phase. 

We conclude that the line profile of OH in RCas is 
produced by the combination of one or two absorption 
lines, one outflow and one infall component, which means 
the occurrence of line doubling at some but not all phases. 
An emission line is running from the blue to the red end of 
the red shifted line. The absorption lines vary in strength 
and slightly in velocity during the cycle and from cycle-to- 
cycle. The emission component can become quite strong, 
so that it becomes visible above the continuum level. A 
weak emission component with an infall velocity seems to 
be present in several spectra. The line profiles of the other 
three miras of our sample can be interpreted in the same 
way. In this way we can also explain the unexpected line 
doubling in x Cyg on JD 2446006 as the passage of an 
emission component through the line core. 

4.2. SiO & HCl 

Representative line profiles for SiO and HCl in the RCas, 
oCet and RLeo spectra are illustrated in Figs.[I], |]and [|, 
respectively. Strong emission components of SiO and HCl 
are seen in several spectra of RLeo and o Cet, and also in 
the spectra of xCyg. In RCas emission components are 
clearly weaker. 

The line profiles of these two molecules show simi- 
lar and complex behavior with both emission and ab- 
sorption. As examples the velocity curves of R Cas and 
X Cyg are discussed (Fig. ||) . One absorption component 
of the SiO and HCl lines in x Cyg at a velocity of about 
+10 km s _1 relative to the systemic velocity of the star 
(about +0.5 kms -1 heliocentric) is visible in all SiO and 
HCl spectra of our data set, including spectra obtained 
outside the time range plotted in Fig.pl. In addition there 
is in xCyg a second absorption component at several 
phases shifting its velocity towards that of the first compo- 
nent. The period observed in these data is about 320 days, 
conspicuously less than the photometric period. Note that 
the velocity curves of the two cycles of x Cyg observed are 
not parallel. We will discuss this difference in more detail 
in the next section. 

R Cas has quite different velocity variations compared 
to the other miras studied. From Figure El we see that 
from JD 2445987 to JD 2446192 the outflow component 
of the SiO lines is dominant in the spectrum, while from 
JD 2446454 to JD 2446570 the infall component is clearly 
stronger. Yet the absorption components have similar 
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Fig. 8. Heliocentric velocities measured from absorption 
components of SiO (filled circles), HCI (open triangles), 
water absorption (stars) and emission (crosses) lines in 
xCyg (upper panel) and RCas (lower panel). The dotted 
line indicates the center-of-mass velocity. 



depths at similar phases. In the spectra obtained in 1986 
( JD 2446430 through JD 2446800) the outflow component 
may be filled by an emission line which is becoming visible 
in the JD 2446570 spectrum. 

A strong emission component is clearly visible in sev- 
eral spectra of R Leo (see Fig.||), oCet and xCyg (Fig.|J). 
The RLeo data set is an example of spectra where the 
line profiles can be interpreted as a combination of an ab- 
sorption line and an emission line, with the absorption 
and emission having different phase dependent velocity 
behavior. On the other hand, in R Cas an additional weak 
emission component at the systemic velocity is seen as an 
obvious asymmetry in the absorption line. In the R Cas 
spectra near maximum light, and especially the spectrum 
of JD 2446192, the strong absorption component of SiO 
and HCI matches the velocity of the strong (for RCas) 
OH emission line. 

As the variability cycle of the 4 /im molecular features 
does not always follow the photometric period it is very 
difficult to compare the line profiles of different miras at 
similar visual phases directly. However, basic characteris- 
tics, namely outflow and infall components, intensity vari- 
ations, and the occurrence of emission at some phases, are 
found in all miras of our sample. The similarity between 
the SiO line profiles and the HCI line profiles is a common 
feature of all stars. In general, line profile and variability is 



similar in RLeo, o Cet and x Cyg- R Cas shows a different, 
more complex behavior. 

- 4.3. H 2 

The H2O line profiles do not follow the variations of either 
the OH or the SiO and HCI line profiles. This is illustrated 
by comparison of H2O line profiles with those of other 
molecules in Figs.p], and EL The dissimilarity of pro- 
files indicates that, at least at some phases, the H2O lines 
probe a very different atmospheric zone than the lines of 
the other molecules seen in the 4/xm spectra. In addition, 
we note that the H2 O lines in x Cyg are weak in agreement 
with the S-type spectral class of this star (Fig.g). Emis- 
sion components, which can be seen in other molecular 
lines, are even more prominent in H2O line profiles. 

Velocities derived from absorption and emission com- 
ponents of the H 2 line profiles of RCas and xCyg are 
plotted in Fig.|| together with the results from SiO and 
HCI. The x Cyg H 2 absorption lines show a velocity vari- 
ation similar to the other three molecules monitored in this 
stars spectrum. The velocity curve is slightly offset from 
the SiO and HCI variations, mimicing the velocities of the 
OH and CO lines but with smaller amplitude. In the one 
case where the CO and H2O velocity significantly differ 
the cause is probably a blend of the emission and absorp- 
tion line profiles. The emission component stays almost 
constant in velocity close to the center of mass velocity 
throughout the entire light cycle. 

R Cas velocities (lower panel of Fig.||) again have quite 
different phase dependent behavior than those of x Cyg- 
One H2O absorption component shows a S-shaped discon- 
tinuous velocity curve apparently spanning the entire time 
of observations, i.e. at least two light cycles! Additional 
absorption components resemble the velocity behavior of 
SiO and HCI around JD 2446200. The emission compo- 
nent seen in x Cyg close to the center-of-mass velocity is 
not found. Emission is instead seen at an infall velocity 
first and switches to an outflow velocity shortly after the 
visual light maximum at JD2446200. The overall behavior 
of the 4 /im spectral lines in R Cas is outstanding and will 
be discussed below. 



5. Mira Dynamics 

Monitoring of the velocity variations of spectral lines pro- 
vides a powerful tool for studying the physical charac- 
teristics and dynamics of the stellar atmosphere. In this 
section we summarize observed velocity variations for the 
two miras with the most complete time series, x Cyg and 
R Cas. On this basis we discuss the dynamics of the outer 
layers of an AGB star's atmosphere. 

There are two basic features of the 4 /xm spectra that 
need to be explained by any model: the occurrence of an 
emission component and the deviation of the variability 
from the photometric period. However, we stress that even 
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a qualitative description of these phenomena has to be 
limited to some basic statements due to the complex, not 
fully understood structure of AGB atmospheres. It can be 
difficult to discern observationally the movement of a layer 
and a change in the optical depth since they can produce 
interfering results. The emission component may have its 
origin in the atmospheric dynamics or in a geometrical 
effect. We will discuss the absorption components first and 
come back to the emission lines at the end of this section. 



5-1- xCyg 

The x Cyg spectra described above provide the most ex- 
tensive and well sampled time series of 4 /mi spectra of 
any star in our sample. \ Cyg also has good time series 
spectra in the H-band and K-band which cover the first 
and second overtone lines of CO. We will summarize the 
velocity data for this star as an example of the overall 
relation between the dynamics of atmospheric layers as 
probed by different molecular features in the infrared. 

The second overtone CO lines at 1.6 /zm are thought to 
provide information about the kinematics deep inside the 
atmosphere, i.e. close to the pulsation driving zone (Hin- 
kle 1978 ). The corresponding velocity curve is S-shaped 
with line doubling occurring around maximum in all mi- 
ras. The line doubling is seen only in miras, not in the 
lower amplitude semiregular variables. In miras the center- 
of-mass velocity is crossed by the velocity curve typically 
around phase 0.4. At that phase no line doubling occurs in 
these stars. A similar behavior is shown by the high exci- 
tation first overtone CO lines found in the 2.3 /mi region. 
For x Cyg an extensive time series of these lines has been 
published already (Hinkle et al. 1982). For the present in- 
vestigation a new set of 1.6 /im spectra was taken in order 
to compare the 4 /im results with simultaneous data from 
CO lines. The 1.6 /im spectra were obtained with the same 
instrumental set up as the earlier measurements published 
in Hinkle et al. ( |l982| ). 

Figure y shows the velocities of the different molecu- 
lar species versus time for the mira x Cyg. Visual phase 
relative to the first maximum observed (phase 0.0) has 
been added as label to t he to p axis. As demonstrated 
previously (Hinkle et al. 1982[ ) the CO lines vary with 
the photometric period. Line doubling is visible around 
maximum. There is no difference between the CO velocity 
curve presented here and those obtained earlier by Hinkle 
et al. fll982| ). The 4 /mi OH lines in xCyg follow the be- 
havior of the CO lines. In general the H2O absorption line 
velocities for x Cyg follow the CO velocity curve. 

The velocity curve defined by the SiO lines is clearly 
different from the CO line curve. As mentioned above SiO 
has a more or less constant velocity component at an infall 
velocity of about lOkms -1 . There is no feature in the CO 
or OH profiles that could be correlated with the first ab- 
sorption component of SiO at constant velocity. A smaller 



velocity amplitude is seen in SiO (and also in HC1, see 
Fig. |) lines. 

SiO also has a second absorption component that 
shows strong variability in velocity. But this variation is 
not always quite with the stellar pulsation. During the 
first cycle observed for x Cyg the line forming region of 
SiO decelerates less than the CO and OH line forming 
region, assuming that the lines of a certain species origi- 
nate in the same layer at all phases. During this cycle the 
braking effect is smaller for this layer. In the second light 
cycle SiO has a velocity component that follows the CO 
velocity curve. 

A possible explanation is that mass loss is not contin- 
uous, with material ejected in shells but not during every 
light cycle (e.g. Winters et al. 2000). Several observations 
in the radio range and of post- AGB objects support this 
description of the mass loss. If so, material would be aggre- 
gated at some distance from the star before it is ejected. A 
mass shell reaching this distance during its pulsation will 
be then confronted with more or less material, i.e. with 
an environment with varying density from cycle-to-cycle. 
If there is less material, the layer or a moving wave can 
go faster and further away from the star and in this way 
again would feed the mass loss. The problem with the time 
between phases of similar velocity being shorter than the 
photometric period can in that way be explained. The lay- 
ers traced by the SiO (and also by HC1) are not varying 
on a shorter time scale but with a different acceleration 
during some cycles. After several cycles we can again have 
similar line profiles as observed in o Cet and x Cyg- 

5.2. RCas 

The mira R Cas shows a behavior remarkably different 
from x Cyg, especially corresponding the velocities of OH 
and SiO lines. No simultaneous CO observations are avail- 
able f or th is star. Comparison with older data (Hinkle 
et al. 1984) shows that the velocity curve derived from 



the OH lines has a flatter slope than the CO velocity 
curve. The maximum expansion velocity observed in OH 
is smaller than the one observed in the CO lines. 

The H2O lines in RCas are the only ones found that 
indicate the presence of (periodic) variations in the outer 
layers of the atmosphere with a time scale longer than 
the period indicated by the visual light change. In Fig.[lC| 
we show that the different H2O absorption components 
can be grouped to two nicely shaped velocity curves by 
two hand drawn fits through the data points. The curves 
are numbered '1' and '2', respectively. Curve '1' as noted 
above covers a time scale exceeding the photometric pe- 
riod, while the variation represented by curve '2' happens 
on a time scale similar to the main period. Curve '2' also 
shows the same relation between velocity and light change 
found e.g. for high excitation CO lines, i.e. minimum ve- 
locity at light maximum. However, no line doubling of the 
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Fig. 9. Velocities of different molec- 
ular lines versus time and visual 
phase for \ Cyg- The visual phase 
given at the top axis is relative 
to the visual light maximum at 
JD 2445838. Crosses mark velocities 
from CO lines, triangles SiO veloci- 
ties and open circles OH velocities, 
respectively. 



same type as in the CO lines is observed for curve '2', as 
any second absorption component fits on curve T'. 

Fig.[10| shows the velocities of the OH lines as well. 
It can be seen that they form a velocity curve of similar 
shape as the H2O lines but with an offset. The first part 
of the OH velocity curve shows also the same velocity 
range as the H2O lines while the second velocity curve 
found from the OH lines has a larger range in velocity. 
This may indicate that our time series does not cover the 
whole length of this long period change, so that the time 
scale may be even larger than two times the main period. 

The line doubling of the H 2 lines suggests that during 
some phases we see line forming layers at two completely 
different locations in the stellar atmosphere. The similar- 
ity in the behavior of curve '2' and the high excitation CO 
lines implies that curve '2' represents a layer close to the 
pulsation driving zone. It is interesting to note that we 
did not find any absorption components that would mark 
this layer around JD2446000. 

Based on these results we propose the following sce- 
nario occurring in RCas: Around JD2445800, i.e. at the 
light maximum preceding our time series, a cycle of ve- 
locity variations started that lasted much longer than a 
typical cycle. A similar model as for the SiO line form- 
ing layer in \ Cyg may therefore be assumed, i.e. that the 
variation marked as curve T' in Fig. nffl lasted longer than 
one cycle as material may have been moved further away 
from the star and needed a longer time to return to its 
original position. 

Note that our data do not give any indication whether 
this behavior is occurring with some periodicity or not. 
The CO lines, observed at a different time, did not show 
such strange behavior. In this context it is interesting to 
note that a very bright maximum is preceding this phase of 
strange velocity behavior. Such bright maxima occur from 
time to time with no obvious regularity in R Cas (Mattei 
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Fig. 10. H2O and OH velocities of RCas versus time. 
Filled and open circles mark the two H2O absorption 
components, crosses indicate the velocities of water emis- 
sion lines, and open triangles denote OH absorption lines. 
Two velocity curves, marked '1' and '2', have been drawn 
through the H2O data points. The horizontal line indicates 
the center of mass velocity. See text for more details. 



|2000 ). It seems quite probable that these two phenomena 
are related to each other. At the next light maximum again 
a velocity curve starts, this time its period is following the 
photometric changes. Therefore, only in some cycles such 
strange line profiles may occur. 

However, as the behavior of the SiO lines in RCas 
shows, this scenario cannot explain the dynamics of the 
outer atmosphere completely. The SiO lines obviously vary 
on a time scale longer than the photometric period, but 
from our data there is no indication for any periodic vari- 
ation. 
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5.3. Emission lines 

Emission can in principle be produced by a temperature 
inversion, as it would be provided by a shock front or a 
hot chromospheric layer, or by a geometric effect in an 
extended spherical object (e.g. Hinkle et al. 200C). Both 
mechanisms are possible in miras. The occurrence of shock 
fronts is well known not only from optical spectra but 
also from the velocity curves derived from second overtone 
CO lines in the near infrared (e.g. Hinkle et al. 1982). In 
the optical range emission of hydrogen and several metal 
lines is triggered by the passage of shock fronts. Chromo- 
spheric layers have been proposed from model calculations 
by Willson & Bowen (1986), but no final observational 



proof has been given for their existence in oxygen rich 
AGB stars. Finally, atmospheres of pulsating AGB stars 
are definitely widely extended and can show a shell like 



structure in the outer parts (e.g. Hofner 1999). 

We note that the emission component of some 
molecules changes velocity during the time series but other 
emission components have a near constant velocity. This 
suggests two separate places for the emission. In R Cas we 
observe an emission component in the OH lines changing 
from a high outflow velocity to the center-of-mass veloc- 
ity. If we attribute the emission to a certain moving layer, 
we would see a braking of the movement until it comes to 
a stop relative to the center of the star. Such a behavior 
of the emission has also been seen in emission lines of the 



Balmer series in miras (e.g. Merrill 1946 ). However, the 
emission seen in the H2O lines, most notably in xCyg, 
is constant in velocity. We propose that this emission is 
produced in the outer layers of the extended atmosphere. 
Tsuji ( 200C( ) makes similar arguments for emission seen 
from water in the inner circumstellar shell. The extended 
atmosphere may also explain emission seen in molecules 
other than H2O in SRV's where there is no indication for 
a shock front from second overtone CO lines or hydrogen 
lines (see below). 



6. SRVs 

Semiregular variables are generally believed to have a 
lower pulsation amplitude than the miras and as a result 
SRV spectra are expected to be less affected by dynamic 
phenomena. To study the influence of stellar dynamics on 
the 4/im spectra of SRVs, we examine in this section rep- 
resentative SR spectra from different period and ampli- 
tude regimes. Short period and low amplitude SRVs typi- 
cally have small mass loss rates and are therefore difficult 
to detect in thermal radio CO lines (e.g. Kerschbaum & 



Olofsson 1999). As a result the most reliable information 



on the center-of-mass velocity for long period variables is 
not available for these stars. The discussion of velocities 
relative to the center of mass is therefore limited to the 
subset of stars where an accurate systemic velocity has 
been measured. 



The SRV observations reported here cover only a single 
epoch. Due to the lack of reliable photometric data a rough 
determination of the actual visual phase of the variables 
at the time of the observation was possible only for three 
SRVs of our sample. As a consequence we will focus on 
WCyg and RUCyg, for which both velocity and phase 
determination are quite reliable. 



Figs. 11 and |12| show the SiO and OH line profile of the 
semiregular variables W Cyg and RU Cyg, respectively. 
The SiO lines in W Cyg show an absorption component 
on each side of the systemic velocity with the outflow 
component dominating. Outflow occurs with a velocity of 



6kme 



In contrast the second overtone CO lines show 



on ly out flow velocities over the whole light cycle (Hinkle et 
al. 1997). One of the CO second overtone spectra has been 
obtained only 2 days before the 4 (im observation and gives 
a heliocentric velocity of — 20.0kms _1 , i.e. an outflow ve- 
locity of about 6kms _1 the same as found for the SiO 
lines. Compared with miras at similar phases (e.g. RLeo, 
line profile at the bottom of Fig. 0) the SiO lines in W Cyg 
are weaker. It should be remarked that the SiO lines in 
WCyg are the weakest among all the SR spectra in our 
sample. Comparison with miras is of course difficult due 
to nonperiodicities. The average OH line profile of WCyg 
(bottom panel of Fig. 0) is broad and slightly asymmet- 
ric towards the outflow part. As for miras, one possibility 
is that the doubling of the absorption is produced by an 
emission component filling in the absorption line, although 
no indication for a shock front is given in the second over- 
tone CO line velocity curve for W Cyg. However, the emis- 
sion is strongly supported by the H2O profile which shows 
emission above the continuum with weak absorption on 
each side. The emission is at — 20kms~ 1 heliocentric. 

The SiO and OH profiles in RU Cyg are quite different 
from those of WCyg. A double structure is visible in the 
OH lines but not in the SiO lines. The spectrum indicates 
only outflow of matter. As for W Cyg a spectrum with 
second overtone lines of CO has been obtained almost si- 



multaneously (Hinkle et al. [19971) . The CO lines indicate 
an outflow, too, with a similar velocity. The profile could 
of course result from an emission/ absorption blend. The 
1/3 lines of water can not be detected in RU Cyg. 

The other five SRVs show slightly asymmetric SiO lines 
with the positive velocity side steeper than the negative 
velocity side. None of these objects shows any indication 
of line doubling or emission components in the line core. 
SiO lines are more narrow in gHer, RLyr and RZAri, 
i.e. in the stars with periods below 100 days, than in the 
other SRVs. The OH lines look similar to the SiO lines in 
these three SRVs. For gHer a reliable systemic velocity is 



available (Kerschbaum & Olofsson 1999). Neither the SiO 



lines nor the OH lines show a significant shift relative to 
this systemic velocity. HC1 lines have been detected in all 
SRVs except for W Cyg. They are much weaker in these 
stars than in miras. 
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Fig. 11. SiO (upper panel) and OH (lower panel) line pro- 
files of W Cyg. The dotted line marks the systemic veloc- 
ity. 
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Fig. 12. SiO (upper panel) and OH (lower panel) line pro- 
files of RU Cyg. The dotted line marks the systemic ve- 
locity. 



H2O is detectable in three of the SR stars, RX Boo, 
g Her, and SW Vir. In all three it is seen in emission. In 
SW Vir the H2O line has an inverse P Cyg shape with 
absorption at — 18.4 km s -1 and emission at — 11.3kms _1 



heliocentric. The H2O velocity for RX Boo is — 11.3kms _1 
and +6.8 km s" 1 for gHer (heliocentric). Comparison with 
center of mass velocities from radio observations shows 
that the water emission occurs always close to but not 
necessarily exactly at the systemic velocity of the star. The 
largest difference was found for gHer (CO radio velocity: 
+1.3kms" 1 ). 



7. Conclusions 

1.1. Explaining low resolution spectra 

Considerable variations in the strengths of SiO bands dur- 
ing a light cycle have been observed in low resolution 
spectra of AGB stars. In our high resolution observations 
the lines are convincingly present at all phases, showing 
a presence of SiO in the atmosphere throughout the light 
cycle. Therefore it is obvious that the simplest propos- 
als, e.g. the dissociation of a large fraction of molecules 
near maximum light, do not explain the spectrum. This 
has been suggested to explain the extreme weakening of 
the SiO bands near light maximum observed in low res- 



olution spectra (e.g. Rinsland & Wing 1982). The iden- 
tification of emission components in these lines suggests 
the blending and cancellation of absorption and emission 



as the explanation for the weakening (Bessell et al. 1996, 
Aringer et al. 1999t). It is interesting to note that Aringer 



et al. (200C) found a significant emission in 4/xm range 



of ISO SWS spectra of RCas, which they attributed to 
SiO, while we found only a rather inconspicuous emission 
component in our high resolution spectra. This discrep- 
ancy suggests strong cycle-to-cycle variations which war- 
rant further investigation. 

7. 2. Dynamics of miras 

A surprise in the 4 /an spectra is the large degree of com- 
bination of velocity variations that arise close to the pul- 
sation driving zone with grossly non-periodic variations 
arising in extended atmospheric layers. In general the con- 
tinuous opacity in the 4/Lim region is significantly higher 
than in the 1.6 /xm range. As a consequence there is a ten- 
dency for lines to originate further out. Therefore many 
molecular features seen at 4/xm monitor stellar variabil- 
ity in the outer parts of the atmosphere, especially many 
SiO and H2O lines. A similar increase in the continuous 
opacity toward the blue is well documented and seen in 
the appearance of the visual spectrum. The line profile 
variations seen at 4/um appear similar to those found for 
Ki and Fei lines around 8000 A by Gillet et al. ( |1985| ). 
Part of the periodic and non-periodic variations seen in 
the 4/zm region could also be due to a deviation from 
spherical symmetry. 

The four stars in our sample show both strong sim- 
ilarities and differences in the 4^m lines. Especially for 
R Cas we find a quite different behavior for most of the 
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lines. We conclude that, while the lower regions of the at- 
mosphere where pulsation originates are very similar in 
all the miras, the outer layers of the mira atmosphere are 
much more distinctive. This could be the result of evolu- 
tionary state along the AGB as might be suggested by the 
range of periods. Other possibilities are differences in the 
fundamental properties like mass, metallicity, or C/O. Or 
the behavior is unique for each star resulting from spa- 
tial inhomgcneitics in the atmosphere. Our sample is too 
small to differentiate between these possibilities. However 
we note that R Leo and o Cet have a similar period and 
similar C/O (based on the water and CO line strengths) 
and have more similarity in their 4 /im spectra than with 
the other two miras. x Cyg, although having a similar pe- 
riod to RCas, is an S-type star so differences between the 
line profiles in these stars are expected to be produced 
by the different atmospheric structure which is related to 
chemical composition. 

Fig.0 shows, that the lines can be a tracer of the 
pulsation, like those in the 1.6 /mi spectrum, but that 
in some cycles irregular behavior occurs. Periodicities in 
the outer layers different from the photometric period 
have been suggested from dynami cal m odel atmospheres 
(se e e.g. Fig. 1 in Hofner & Dorfi 1997 , and Hofmann et 



al. 1995). The observational findings can be qualitatively 



understood in the framework of these dynamical models 
of AGB star atmospheres. During some cycles, which seem 
to be marked by a brighter visual maximum, part of the 
outward moving layers decelerates less than in other cycles 
due to a smaller density in the outer layers. This density 
variation may be related to a non continuous mass loss. 
Such deviations are not monitored by lines formed deeper 
inside the atmosphere. It is therefore an effect only occur- 
ing in the outer parts of the star. This means that this 
behavior is not related to a different strength of the pul- 
sation. 

In this co ntext it is interesting to note, that Rinsland & 
Wing ( 1982 ) made two low resolution 4 /mi spectra of each 
of the miras RCyg, x Cyg and RCas separated by almost 
exactly one light cycle. In two miras, R Cas and x Cyg, 
they detected a clear difference in the overall strength of 
the SiO (2,0) and (3,1) bands between the two observa- 
tions made at similar phase. This is in good agreement 
with the findings discussed above. 

7.3. SRVs 

Small amplitude variables like the SRVs of our sample are 
thought to be more compact than miras. Variations of CO 



lines at 1.6 /mi (Hinklc et al. 1997, Lebzelter 1999) indicate 
that SRVs vary with a much smaller velocity amplitude 
and that shock fronts are very weak or do not occur at all. 
In agreement with this finding, the SRVs of our sample do 
not show such extreme emission components as found in 
miras. The line profiles of most SRVs show a similar shape 
and strength. This indicates that their outer atmospheric 



layers show a similar structure. Due to the lack of time se- 
ries for these stars no study on their individual variability 
is possible. 

As our main result for SRVs we note that both the 
occurence of emission components as well as the asymme- 
tries and the shifts relative to the center of mass velocity 
clearly indicate that velocity fields have to be included into 
a model atmosphere to reproduce the results properly. Hy- 
drostatic models cannot be used for high resolution spec- 
tra of these stars. However, the fact that the effects of the 
stellar pulsation are still rather small in SRVs compared 
to miras makes these stars excellent starting points for 
adapting dynamical models to the observations. 



8. Outlook 

The aim of this paper was to present and to discuss qual- 
itatively the observational results of our high resolution 
spectroscopic monitoring of AGB variables. In a forth- 
coming paper we will present synthetic spectra based on 
dynamic model atmospheres in order to model the obser- 
vations presented here. 

However, an obvious explanation of the non-regular 
variations is inhomogeneities in the upper atmosphere. 
These inhomogeneities could be in the form of filaments 
or clouds in the extended stellar atmosphere. Unique fila- 
ment structure would explain non-cyclic behavior and the 
strong spectral identity of individual stars. There is con- 
siderable evidence from maser observations of such con- 
densations in the outer atmosphere. The SiO maser, which 
is probably formed in the same atmospheric region as some 
of the 4 /mi lines, is clearly formed in condensations (Di- 
amond et al. 



1994). Large scale inhomogeneities are indi- 



cated by maps and images of AGB circumstcllar shells. A 
further indication for inhomogeneities in the atmospheres 
of AGB stars come from polarimetry measurements. As 
shown from model calculation (Harrington 1969) obser- 
vational results can be understood if it is assumed that 
temperature variations occur over the stellar surface. This 
model has been applied succ essful ly to observations of o 
Cet by Tomaszewski et al. ( |1980| ). Such filament struc- 
tures currently exceed foreseeable modeling capabilities, 
but it might be possible to include their effects in future 
dynamical model atmospheres. 
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